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Abstract 
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extracted from the frequency dependence of the coupling coefficient for the far-infrared absorption 
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dynamics model of AS2S3 glass qualitatively supports these findings. 
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I. INTRODUCTION 



Disordered solids, such as glasses, are characterized by a lack of structural order, which 
leads to disorder in many physical quantities. For example, atomic charges in glasses fluc- 
tuate in space. The origin and properties of such fluctuations are of considerable interest in 
the field of disordered systems. One such intriguingquestion concerns the scale of charge 
fluctuations and possible charge ordering in glasses [1, . Another intriguing ques- 

tion is related to the possibility of obtaining information about the charge distribution from 
experiment. It has been recently discussed how terahertz absorption spectroscopy can be 
used for this purpose [6j. 

In this paper, we study the atomic charge distributions in glassy arsenic sulfide, As 2 S 3 , 
by means of first principle tight-binding molecular dynamics simulations and compare the 
characteristics of these distributions with those obtained from experimentally available far- 
infrared (FIR) absorption coefficient measurements. 

II. THEORETICAL BACKGROUND 

First, we recall how the absorption coefficient in the FIR region is related to the atomic 

n 

charge distribution |6J. The expression for the linear absorption coefficient of photons, a(co>), 
caused by the interaction with harmonic atomic vibrations in solids is given by the following 
expression (see e.g. jlj]), 

9(f*)) , (1) 
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where and Z iya p are the mass and the dynamical charge tensor for atom i (i = 1, . . . , N, 
with N being the number of atoms in a solid of volume V; a and (3 run over Cartesian 
coordinates), eip(u) is the /3-component of the vibrational eigenvector of frequency uj cor- 
responding to atom i, stands for the high-frequency dielectric constant, n = N/V is 
the atomic concentration and angular brackets denote configurational averaging. As follows 
from Eq. (TjQ), the IR absorption coefficient is proportional to the self- averaging value of the 
vibrational density of states (VDOS), 
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where j runs over all the eigenfrequencies, i.e. 

a(u) = (C(u)) g(u) . (3) 

The coefficient of proportionality, (C(u)), is called the coupling coefficient between IR pho- 
tons and atomic vibrations for linear light absorption. The VDOS behaves universally 
(g(w) oc uj 2 according to the Debye law) in the FIR regime and thus all interesting and 
possibly universal features can be attributed to the frequency dependence of the coupling 
coefficient. 

The absorption coefficient given by Eq. (CQ) is temperature independent and some 
variations of a(u) with temperature found experimentally for microwave frequencies {zj], 
uj/2itc < 1 cm -1 , which may possibly be attributed to excitations of two- level systems 
and/or to highly anharmonic atomic modes, are not considered here. The dynamical charge 
tensors are crucial for a correct description of the peak positions and their relative intensities 
in the bulk of the vibrational band (above the FIR region) S| but are not so significant 
in the FIR regime and thus a simpler (rigid-ion) model can be used, so that 
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C(u) = C„ 
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(4) 



with qi being the temporally fixed but spatially fluctuating atomic charges and 

C = 27r 2 n/cyfeZ ■ (5) 

In the well-studied case of ordered systems, where the charges do not fluctuate, and the 
eigenmodes are phonons, the coupling coefficient is non-zero only for optic modes at the 
centre of the Brillouin zone. In disordered systems, structural disorder leads to charge 
transfer between atoms, i.e. to disorder in atomic charges gj, and to intrinsic disorder in the 
components of the eigenvectors which lose their translational invariance. These two related 
sources of disorder, encoded in Eq. are responsible for the peculiar behaviour of C(u) 
in amorphous systems, which actually has a universal functional form in the FIR regime |6( , 

(C(u)) ~ A + Bu 2 , (6) 

where A and B are material-specific constants. 

In order to see this, we use two facts known about the structure of the eigenmodes in the 
FIR regime and about the distribution of atomic charges. First, the disordered eigenmodes 



in the FIR regime resemble plane waves 10|, lU( characterized by pseudo-wavevectors k and 



exhibiting pseudo-dispersion, uiAk) = Cjk (with Cj being the sound velocity for branch j) 
and can be well approximated [6| by a plane wave characterized by wavevector k and unit 
polarization vector p^, 

with fn = N~ 1 £\ mi and being the position vector of atom i. 

The second useful fact concerns the distribution of charges in disordered systems. It 
has been found in simulations that the charges in the models of some amor- 

phous materials preserve approximately charge neutrality within certain structural units. 
7or example, the SiC>4 structural units in vitreous silica are approximately electro- neutral 
l| meaning that the positive charge on an Si atom is approximately equal in magnitude 
to half of the sum of the charges on the four nearest oxygen atoms. The values of silicon 
and oxygen charges vary strongly between structural units, depending on local structural 
characteristics such as the Si-O-Si bond angle. Moreover, the electro-neutrality within the 
structural units is maintained only approximately (see below) and there is always a stochas- 
tic component in the charge distribution due to intermediate and long-range fluctuations in 
the structure. These observations allow the values of atomic charges to be split into two 
components, qi = qu + q 2i , with (^({rj}) representing uncorrelated charge components, and 
the random charges q 2i satisfying local charge neutrality. 

The values of qu({ri}) depend on the atomic coordinates {r^} in a complicated fashion 
so that we can approximately assume the absence of correlations between qu on different 
atoms, i.e. 

(quqij) ~ (qu)(qij) ^ cr?^ , (8) 

where the variance a\ { can vary for different type of atoms, or for the same atoms but 
e.g. that are abnormally coordinated. Similarly, we assume no correlations between q u and 
atomic position vectors, so that 

( gi ,e ik ^)~( gi ,)( e ik ^)~0. (9) 

The random charges q 2 % obeying local charge neutrality can be imagined as resulting 
from charge transfers between nearest-neighbour atoms, i.e. q 2 % = where j runs 

through all the nearest-neighbours of atom i and Aqji (= —Aqji) is the charge transfer from 
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the originally neutral atom j to the originally neutral atom i. In heteropolar crystals, the 
values of Aqji are not random and finite. In disordered systems, the values of Aqji are 
distributed around mean value (s) which do not necessary coincide with those for crystalline 
counterparts (see e.g. Such fluctuations and deviations of means in Aqji are due to 

distortions in local structural units, e.g. in bond angles and bond lengths. We also assume 
that there are no correlations between randomly fluctuating charges qu and local charge 
transfers, 

(Qu92j) ^ (Qu)(Q2j) - ■ (10) 

The configurationally averaged coupling coefficient can be recast in terms of correlated 
and uncorrelated charges in the following manner, 



Cn 



E^ elk ' Ri + E^ ik - Ri 
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where S n = ^ fce ±,R ' (n = 1,2). In the absence of correlations between qu and q 2 j (see 
Eq. (JTOl ). the above formula for C(u) reduces to 

C(u) = N- 1 ((\S l \ 2 ) + (\S 2 \ 2 )) . (12) 

The first component in Eq. (TT2"]) can be further simplified as, 

iV- 1 (l^l 2 ) * iV" 1 5> 4gy ) <e ik -^-)) ^ n-i 4 = ^ , (13) 

ij i 

where we have used Eqs. OH])-©- In the case of a two-component system containing Ni and 
A^ 2 atoms of different types, a\ = (Ni/N)^^ + (Ni/N)af 2 - Therefore, the first contribution 
in the coupling coefficient is frequency independent and depends only on the variance of 
uncorrelated charge distributions. 

The second component in the coupling coefficient, iV _1 (|5*2| 2 ), which is due to random 
and locally correlated charge fluctuations, does not contain the frequency-independent part 
and, in fact, is proportional to u 2 . This can be demonstrated using the bond representation 
for 5*2, 

S2=[Ys A %' ( eik rj - ) = ( A % eik ' r(lJ) (e ik ' r " /2 - e ik ^ /2 ) \ , (14) 

with Vij = Yj — Yi and Xuj\ = (rj + Yj)/2, where the sum is taken over all the bonds (ij) in 
the system. In the absence of the plane wave (k=0), this sum equals zero, thus reflecting 



global charge neutrality of locally neutral units (the contribution from each bond is exactly 
zero due to the local charge neutrality). In the FIR regime, k • <C 1 and Eq. (TH|) can be 

recast as 

S 2 ~kJ2 A <lve ik ' T(ij) (™-rij) , (15) 

with k = kh. Consequently, the contribution from the correlated charges to the coupling 
coefficient is 

(\S 2 \ 2 )^k 2 / J2 AgyAfty^W)-^)^.^)^-^)^ , (16) 
\(ij)(i'f) I 

so that (\S 2 \ ) oc k 2 oc uj 2 (implying linear dispersion in the FIR regime). The double sum 

in Eq. ffl6|) depends on precise structural details of the material but does not depend on k 

in the FIR range (oc Const + 0{k 2 )) and thus the (^-dependence of the second contribution 

in the coupling coefficient is a general feature of the FIR absorption in disordered solids. 

Therefore, we have demonstrated that uncorrelated atomic charges result in the 

frequency-independent part of the coupling coefficient, while locally correlated charges, 

maintaining the charge neutrality within local structural units, give rise to the quadratic 

frequency dependence of the coupling coefficient. 



III. RESULTS 



As follows from the previous section, the coupling coefficient for the FIR absorption has 
a universal frequency dependence (see Eq.pl)), containing a frequency- independent part and 
a frequency- dependent contribution (oc uj 2 ). The frequency dependence of the absorption 
coefficient, ct(u), can be measured experimentally using THz time-domain spectroscopy 
B, Q. The VDOS cau also be measured experhueutahy using, e.g. .elastic ueu tr ou 
scattering [131 . 114J ]. Therefore, the frequency dependence of the coupling coefficient for the 
FIR absorption, C{uj) = a(u>)/g(uj), can be found experimentally (see Fig. [TJ). Fitting 
experimental data with theory ([6]) allows the values of constants A and B entering Eq.([6]) to 
be estimated. We have done such a fitting to the experimental data for AS2S3 glass (cf. the 
solid and dashed lines in Fig.[TJ) and found the values of these constants to be A = 1780 cm -2 



and B = 75 cm -1 

The constant B depends on the structural characteristics of glass in quite a complicated 
fashion (see Eq. fflBI) and cannot be used straightforwardly for extracting the relevant charge 



distributions. However, the constant A is directly related to the width of the uncorrelated 
charge distribution, 

1 = ~7y ' I 17 ) 
^0 

where the coefficient Co is given by Eq. ([5]). Using Eq. ( JT7j) . we have estimated the value of 
o 7 ! ~ 0.12 e for a-As2S3. 

In order to verify and support the consistency of the presented model for the frequency 
dependence of the coupling constant, we have created a model of As 2 S 3 glass using de nsity 
functional theory-based tight-binding (DFTB) Q] molecular-dynamics simulation 17]. 
First, we checked the local charge neutrality within AsS^ pyramids, these being typical struc- 
tural units in AS2S3 glass. The results presented in Fig. [2] clearly demonstrate correlations 
between the charge at the central As atom and surrounding S atoms. As expected for such 
a covalent material, the local charge neutrality is not exact and the data points fluctuate 
around the bisector due to contributions from qu. 

Second, we calculated the Mulliken charge distributions for As and S atoms in our DFTB 
model of AS2S3 glass (see Fig. [3]). The charges are distributed approximately normally (see 
the dashed lines in Fig. [3])) around the mean values q As ~ 0.56 and q s ~ —0.36, with the 
standard deviations ai As ~ 0.06 and o"is — 0.05. Therefore, the value of a can be estimated 



as Wi = a/(2/5)ct 2 As + (3/5)cr 2 As ~ 0.054. This value is somewhat less than that estimated 
from the fit of the experimental data for C(u) by Eq. ([6]) (cf. the solid and dashed lines in 
Fig. []]), i.e. a 1 ~ 0.12. Several possible effects could account for such a discrepancy. The 
Mulliken charges are an artificial way of assigning charge values to particular atoms within 
the DFTB scheme - a different population analysis may give different absolute values of the 
charges and thus different widths of the charge distributions. Another possible reason for the 
discrepancy is due to the experimental uncertainty in measuring the absorption coefficient 
at very low frequencies, u < 5cm _1 , caused by parasitic secondary reflections of THz pulses 
of light which, after the Fourier transform, can contribute in this frequency range. 



IV. CONCLUSIONS 



To conclude, it has been demonstrated that the frequency dependence of the coupling 
coefficient for far-infrared absorption can be used to extract characteristics of the atomic 
charge distribution in glasses. Namely, a fit of such an experimentally measured dependence 
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FIG. 1: (Color online) Experimental frequency dependence of the absorption coupling coefficient 



la]). The dot-dashed lines represent the 



in the FIR range for AS2S3 (solid line 6J] and circles 
numerical data obtained from a DFTB molecular-dynamics model. The dashed line shows the fit 
of the experimental data by Eq. [6] with A = 1780 cm -2 and B = 75 cm . 

by the theoretically predicted law, A + Bu 2 , allows the constants A and B to be extracted. 
The value of A is simply proportional to the variance of the uncorrelated charge distribution, 
o\. In the case of AS2S3 glass, we have estimated this value to be eFj ~ 0.12. A similar 
estimate of a\ ~ 0.054 has been obtained from a first-principles molecular-dynamics model 
of the same glass. 
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